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Zigzag graphene nanoribbons patterned on graphane are studied using spin-polarized ab ini- 
tio calculations. We found that the electronic and magnetic properties of the graphene/graphane 
superlattice strongly depends on the degree of hydrogenation at the interfaces between the two 
materials. When both zigzag interfaces are fully hydrogenated, the superlattice behaves like a free- 
standing zigzag graphene nanoribbon, and the magnetic ground state is antiferromagnetic. When 
one of the interfaces is half hydrogenated, the magnetic ground state becomes ferromagnetic, and 
the system is very close to being a half metal with possible spintronics applications whereas the 
magnetic ground state of the superlattice with both interfaces half hydrogenated is again antifer- 
romagnetic. In this last case, both edges of the graphane nanoribbon also contribute to the total 
magnetization of the system. All the spin-polarized ground states are semiconducting, independent 
of the degree of hydrogenation of the interfaces. The ab initio results are supplemented by a simple 
tight-binding analysis that captures the main qualitative features. Our ab initio results show that 
patterned hydrogenation of graphene is a promising way to obtain stable graphene nanoribbons with 
interesting technological applications. 

PACS numbers: 73.22.-f, 73.22.Pr, 71.15.Mb 



I. INTRODUCTION 

Graphene^ has emerged as a prospective material for 
nanoelectronics^— and spintronics applications,^ 7 owing 
to the high mobility of carriers, the low-intrinsic spin- 
orbit interaction, as well as the low-hyperfine interaction 
of the electron spins with the carbon nuclei^ The absence 
of an energy bandgap in graphene can be overcome by 
patterning graphene into nanometer size ribbons, which 
have been realized in recent experimental studies^ - — The 
existence of width-dependent energy bandgaps makes 
graphene nanoribbons (GNRs) a potentially useful struc- 
ture for various semiconducting applications. Ab initio 
calculations show that zigzag GNRs are magnetic^ and 
they become half metallic by breaking the symmetry be- 
tween the left and right edges <2r— This can be achieved 
by using a different chemical termination at each edge^iiii 
or by applying an external electric field. 11 

Transport properties of GNRs depend on its edge 
disorder A At present, controlling the shape and width 
of the GNR is a very challenging experimental task4£ 
There are also experimental evidences, corroborated by 
ab initio calculations, that the zigzag edge is metastable 
in vacuum because it shows a planar reconstruction to 
lower its energyJ^ A new possibility for the fabrication 
of high-quality graphene nanoribbons is by selectively hy- 
drogenating graphene or by carving GNRs on a graphane 
sheets 

First-principles calculations predicted that graphane 
is a wide-gap insulator! 14 : 15 This novel material 
was synthesized recently by exposing graphene to a 
cold hydrogen plasma^ The promising properties of 
graphene/graphane structures have attracted recently 
a large interest ! 13 : 17 " — It was shown that the energy 



bandgaps of both zigzag and armchair graphane nanorib- 
bons increase as the nanoribbons become narrower! 17 : 18 
Hybrid graphene/graphane nanoribbons in vacuum 
where also studied from ab initio calculations and it 
was shown that the energy bandgap is dominated by the 
graphene section of the nanoribbon^ A bandgap open- 
ing in graphene, induced by the patterned absorption of 
atomic hydrogen onto the Moire superlattice positions 
of graphene grown on a substrate, was recently found 
experimentally. 19 Ab initio calculations indicate that the 
observed gap opening is due to the confinement effect in 
the residual bare graphene regions^ 

The edge states of graphene/graphane superlattices 
where recently considered in the tight-binding approx- 
imation. If the interface is oriented along a zigzag di- 
rection, edge states enhance effects related to the spin- 
orbit interaction^ The stability of graphene/graphane 
nanoribbons was also studied from ab initio calcula- 
tions, and it was predicted that the graphene/graphane 
nanoribbons are stable down to the limit of a single car- 
bon chain<2i 

Patterned fluoridation of graphene is another poten- 
tial way for making graphene nanoribbons. Experi- 
mental evidences of graphene fluoridation was recently 
reported! 22 ' 23 Graphene fluoride is a two-dimensional ma- 
terial that exhibits an insulating behavior similar to 
graphane; 22 : 23 which can be also used to create super- 
lattices of graphene nanoribbons. 

In this paper, we carry out spin-polarized ab initio cal- 
culations to study the electronic and magnetic properties 
of superlattices of zigzag graphene/graphane nanorib- 
bons (ZGGNRs). We study three cases: (i) when all 
the interfaces between the graphene/graphane nanorib- 
bons are fully hydrogenated, (ii) when one edge of the 
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nanoribbons is fully hydrogenated and the other edge 
is half-hydrogenated, and (iii) when all the interfaces of 
the superlattice are half hydrogenated (see Fig. [T]). We 
found that the magnetic ground state of the ZGGNRs 
are semiconducting in the above three cases but the elec- 
tronic and magnetic properties strongly depend on the 
degree of hydrogenation of the interfaces. 

The calculations were performed with the QUANTUM 
ESPRESSO package- 4 employing spin-polarized density- 
functional theory and the Perdew-Burke-Ernzerhoff 
exchange-correlation functional^ An ultra-soft pseu- 
dopotential description of the ion-electron interaction^ 
was used together with a plane-wave basis set for the 
electronic wave functions and the charge density, with 
energy cutoffs of 40 Ry and 400 Ry, respectively. The 
electronic Brillouin zone integration was sampled with 
a 16 x 1 x 1 uniform fc-point mesh and a Gaussian 
smearing of 0.01 Ry. The two-dimensional behav- 
ior of the graphene/graphane nanoribbons was simu- 
lated by adding a vacuum region of 12 A above the 
graphene/graphane sheets. All the structures were re- 
laxed using a criteria of forces and stresses on atoms of 
O.OOleV/A and 0.5GPa, respectively. 



II. RESULTS 

One or both of the two carbon atoms in each unit cell 
of the zigzag chain at the graphene/graphane interface 
can be hydrogenated (see Fig. [1]). The degree of hydro- 
genation of the zigzag chains has a very important in- 
fluence on the electronic properties of the ZGGNRs. In 
the following, we will refer to the fully covered hydro- 
gen interface as an a interface and to the half-covered 
interface as a /3 interface, as shown in Fig. [TJ In prin- 
ciple, if we allow the hydrogen deposition to start on a 
certain sublattice of graphene, it would be possible to ex- 
perimentally control the degree of hydrogenation of the 
interfaces. 

Zigzag graphene nanoribbons (ZGNR) are usually clas- 
sified by the number of zigzag chains across the ribbon 
width&ii In order to simplify the analysis, we consid- 
ered superlattices of zigzag graphene/graphane nanorib- 
bons where both, graphene and graphane nanoribbons, 
are six-chains wide. This means that the unit cell of Fig. 
[TJa) has 24 carbon and 12 hydrogen atoms. The degree of 
hydrogenation at the interface, and therefore the number 
of hydrogen atoms, decreases from Fig. QIa)-QTc). 

The interface between graphene and graphane nanorib- 
bons is an interface between sp 2 and sp 3 bonded carbon 
atoms. Carbon atoms in graphene have a sp 2 hybridiza- 
tion and the structure stabilizes as a flat sheet. The 
hybridization changes from sp 2 to sp 3 in graphane due 
to the presence of alternating hydrogen atoms placed on 
both sides of the carbon plane. The sp 3 hybridization 
also forces nearest-neighbor carbon atoms in graphane to 
lie in different planes, and the resulting graphane struc- 
ture is undulating. 15 
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FIG. 1: (Color online) Relaxed zigzag graphene/graphane 
nanoribbons. (a) Both graphene/graphane interfaces are fully 
hydrogenated (two a interfaces), (b) one of the interfaces is 
fully hydrogenated and the other one is half-hydrogenated (a/3 
nanoribbon), and (c) both interfaces are half-hydrogenated 
resulting in nanoribbons. Carbon atoms are represented 
by large (yellow) spheres and hydrogens atoms by small (light 
blue) spheres. We include a side view and a perspective view 
of the relaxed structures for each case. The unit cells of the 
superlattices are shown by dashed-dotted lines in the side 
views. 



Relaxed structures of the ZGGNRs are shown in Fig.[T] 
All structures were fully relaxed according to the forces 
and stresses on the atoms by using spin-polarized cal- 
culations. When both interfaces are fully hydrogenated 
(two a interfaces) one side of the graphene nanorib- 
bon is pushed up and the other side is pushed down by 
the graphane sp 3 -bonded carbon atoms. This induces a 
tilt angle of 164° between adjacent graphene/graphane 
nanoribbons as shown in Fig. [Ha). To satisfy the 
conditions imposed by the sp 3 / sp 2 carbon interface, 



the graphene nanoribbon is tilted with respect to the 
graphane plane but it remains flat. A similar situation 
occurs when both interfaces are half-hydrogenated (two 
/3-type interfaces) as shown in Fig.QJc). In this case, the 
tilt angle between the graphene and graphane nanorib- 
bons is 166.6 degrees. When one interface is fully hy- 
drogenated and the other one is half-hydrogenated (a/3 
superlattice), the graphene nanoribbon is pushed down 
at both edges by the graphane sp 3 -bonded carbon atoms. 
This results in a curvature of the graphene nanoribbons 
as shown in Fig. [IJb). 

It is known that a freestanding a/3 graphene nanorib- 
bon is not stable. We found that by using spatially se- 
lective hydrogenation of graphene we are able to create 
stable a/3 nanoribbons. The stability of those nanorib- 
bons was checked by allowing relaxation of the structure 
in our ab initio calculations. 

We also addressed the question about the feasibility of 
appearance of the zigzag interface by comparing it with 
the armchair interface. We prepared two supercells (one 
with a zigzag and the other one with an armchair arrange- 
ments of atoms) forming graphene/graphane nanorib- 
bons and containing the same number of atoms (108 
atoms), 72C and 36H. The zigzag supercell is equiva- 
lent to three times the unit cell of Fig. HJa) replicated 
along the graphene/graphane nanoribbon. We obtained 
a very small energy difference between these zigzag and 
armchair graphene/graphane nanoribbons, only 0.002 
eV/atom. This value is smaller than our numerical ac- 
curacy, and it is also very small if we compare it with 
the thermal energy at room temperature, 0.025eV. We 
can conclude that the formation energy of zigzag and 
armchair graphene/graphane nanoribbons are equally fa- 
vorable at room temperature. This situation is different 
from the case of freestanding nanoribbons where it was 
found that the zigzag edge is metastable, because it shows 
a planar reconstruction to lower its energy^ This shows 
that graphene/graphane superlatticcs are a more conve- 
nient way for stabilizing zigzag graphene nanoribbons. 

In calculations not shown here, we carried out struc- 
ture relaxations where the volume of the unit cell is 
kept constant and equal to the volume of an equiva- 
lent graphane system and where only the forces between 
atoms are relaxed. This is equivalent to the experimen- 
tal situation where graphene nanoribbons are carved in 
a graphane sheet with fixed borders. In this case, we 
found no tilt angles between the graphene and graphane 
nanoribbons and only a very small curvature for the a/3 
superlattice. The electronic or magnetic properties of the 
system did not change in a significant way if we compare 
with the fully relaxed case. 

In the following sections, we discuss separately the elec- 
tronic and magnetic properties of the ZGGNRs with dif- 
ferent degree of hydrogenation at the interfaces. 





FIG. 2: (Color online) Three-dimensional (3D) graph of the 
spin density (p-f — p±) of the act zigzag graphene/graphane 
nanoribbon. (a) Spin density of the antiferromagnetic ground 
state and (b) spin density of the metastable ferromagnetic 
state. In both cases, the magnetization is almost completely 
localized in the graphene nanoribbons, and the edges of the 
graphane nanoribbons are only slightly magnetized. Regions 
with positive spin density are represented in red (dark gray) 
and regions with a negative spin density are shown in light 
blue (light gray). Isosurfaces corresponds to (pf — pi)—0.02 
e/A 3 , where e is the electron charge. 



A. Two fully hydrogenated interfaces: aa zigzag 
graphene/graphane nanoribbons 

We first analyze the case where both zigzag interfaces 
are fully hydrogenated. The ground state of the sys- 
tem corresponds to an antiferromagnetic (AFM) order- 
ing between the two opposite graphene edges and a ferro- 
magnetic (FM) ordering along the zigzag carbon chains. 
We will refer to this state as an AFM state. The spin 
density (p^ — Pi) of the AFM state is represented in 
Fig. Ufa) while the metastable FM state (ferromagnetic 
ordering both between the two graphene edges and along 
the zigzag carbon chains) is represented in Fig. &h). 

The total magnetic moment and the absolute mag- 
netic moment for the FM state are 0.59/is/cell and 
0.95/xs/cell, respectively (ps is the Bohr magneton). 
The AFM state has an absolute magnetic moment of 
1.09/XB/cell. 

The spin density (p-f— p\) of the aa ZGGNR are mostly 
concentrated in the graphene nanoribbons as shown in 
Figs.JHa) and[UJb). When moving away from the inter- 
face to the center of the graphene nanoribbons, the spin 
density decreases faster in the FM state than in the AFM 
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FIG. 3: (Color online) The spin-resolved band structures of 
the superlattice with all the interfaces fully hydrogenated, 
aa ZGGNR superlattice. (a) Antiferromagnetic semicon- 
ducting ground state, (b) ferromagnetic metallic metastable 
state, and (c) nonmagnetic metallic metastable state, (d) The 
band structure of the AFM state of a hydrogen-passivated six- 
chains wide zigzag graphene nanoribbon. In all figures, the 
Fermi energy is set to zero. The energy differences of the 
ferromagnetic (b) and the nonmagnetic states (c) with re- 
spect to the antiferromagnetic state (a) are Eafm — Efm = 
8meV/edge atom and Eafm ~ Enm = 35meV/edge atom, 
respectively. Spin up and spin down bands coincide in (a), 
(c), and (d). 



state. As a consequence, in the middle of the graphene 
nanoribbon, the spin density is lower in the FM state 
than in the AFM state. This explains the difference of 
0.14^s/cell in the absolute magnetization between both 
magnetic states. 

We also carried out ab initio calculations with a unit 
cell twice larger along the zigzag carbon chains to explore 



TABLE I: Comparison between a six-chains wide freestanding 
zigzag graphene nanoribbon (ZGNR) and a six-chains wide 
aa zigzag graphene/graphane nanoribbon (aa ZGGNR). 
Energies are in millielectron volt per edge atom, bandgaps in 
millielectron volt, and the absolute magnetic moments (i^abs) 
in fiB /cell. The total magnetic moments for the FM state 
are 0.51^ fl /cell and 0.59^ fl /cell for ZGNR and aaZGGNR, 
respectively. 
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whether an antiferromagnetic ordering along the nanorib- 
bons can be the ground state of the system. However, 
we found that ferromagnetic order along the nanoribbon 
shown in Figs. [Ha) and [H[b) gives always lower energy 
states. 

The band structure of the ZGGNRs with the two inter- 
faces fully hydrogenated are shown in Fig.[3J Figure 03a) 
corresponds to the semiconducting AFM ground state. 
The FM metastable state [Fig. [3jb)] shows a metallic 
behavior. As the graphane nanoribbons are not spin po- 
larized there is no significant spin density observed that 
can be associated with the graphane regions. The non- 
magnetic metastable state is a semimetal as shown in 

Fig.ETc). 




FIG. 4: (Color online) Side and perspective views of spin 
density of the ferromagnetic ground state of the a/3 ZGGNRs. 
Around the j3 interface the spin density is larger than in the 
a interface. The /3-edge of the graphane nanoribbon is also 
magnetized, while the spin-state of the a edge is similar to 
Fig. [2] Regions with positive spin density are represented in 
red (dark gray) and regions with a negative spin density are 
shown in light blue (light gray). Isosurfaces corresponds to 
( Pt - Pi )=0.02 e/A 3 . 

The band structure of the AFM state of an equiva- 
lent six-chains wide ZGNR, thus without graphane, is 
represented in Fig. [3jd). We see that, near the Fermi 
level, the band structure of the AFM state of the super- 
lattice [Fig. [3] (a)] is very similar to the band structure 
of the AFM state of a graphene nanoribbon in vacuum 
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[Fig. [3]Jd)]. This similarity is also observed in all the 
other metastable states. Thus, near the Fermi level, the 
band structure of the superlattice is dominated by the 
contribution of the graphene nanoribbons. 

We carefully compared the magnetic behavior of iso- 
lated six-chains wide ZGGNR with graphene nanorib- 
bons embedded in graphane. The comparison between 
their electronic and magnetic properties are summarized 
in Table HI We found very small differences in the energy 
bandgaps, magnetic moments, and energy differences be- 
tween the different magnetic states of graphene nanorib- 
bons in vacuum and patterned in graphane. These results 
show that patterning graphene nanoribbons in graphane 
is a promising way of obtaining zigzag graphene nanorib- 
bons with the same electronic and magnetic properties 
as in vacuum. 

In the four cases represented in Fig. [31 localized edge 
states generate flat bands near the Fermi level. Infor- 
mation about the spatial localization of the edge states 
can be obtained by analyzing the partial (or band de- 
composed) charge density (PCD). This magnitude is also 
called integrated local density of states (ILDOS), which 
describes the density of states but space resolved, and it 
is useful when interpreting the data from STM experi- 
ments. In the following, we will use the term PCD. 

Figure Ufa) shows the spatial localization or PCD of 
the edge states around Ef (Ef ± 0.02eV) for the non- 
magnetic case. If we exclude the states at Ef, and do an 
integration from Ef — 0.09eV to Ef — leV, we see that the 
states become less localized, see Fig.[5jb). An interesting 
feature in Fig. [SJb) is that near each edge of the graphene 
nanoribbon only one of the two carbon sublattices makes 
a large contribution to the states near the Fermi level. 
In fact, we can appreciate from Fig. EJb) that different 
carbon sublattices contribute at both edges while just in 
the middle of the graphene nanoribbon both sublattices 
contribute to the bands. The main contribution of the 
other carbon sublattice is located near the center of the 
Brillouin zone at the T point. These states lie at lower 
energies {Ef — 2eV) and are entangled with the bands 
coming from the graphane nanoribbon. 

All the above mentioned features of the band structure 
and the PCD of aa ZGGNRs can be captured qualita- 
tively with a simple tight-binding model, which results 
are summarized in Figs. O^a) and^d) of the Appendix. 



B. Superlattice with one interface fully 
hydrogenated and the other one half-hydrogenated: 
a/3 ZGGNRs 

When one of the interfaces is half-hydrogenated and 
the other one is fully hydrogenated the spin-polarized 
ground state of the system corresponds to a ferromag- 
netic ordering, both along the graphene/graphane inter- 
face and between the two edges of the graphene nanorib- 
bon. Figure 2] shows the spin density for the FM ground 
state of the a/3 ZGGNR. 




FIG. 5: (Color online) 3D graph of the PCD. (a) PCD of 
the nonmagnetic state of an aa superlattice integrated in the 
energy window Ef ± 0.02eV, (b) the same as before but from 
Ef — 0.09eV to Ef — leV, (c) PCD of the majority spin states 
of an a/3 superlattice in the FM state integrated in Ef ± 
0.7eV, (d) the same as (c) but for the minority spin states, 
and (e) PCD of the AFM state of a (3(3 superlattice from Ef 
to Ef — 0.16eV. Red regions corresponds to majority spins (c), 
blue to minority spin (d), and violet means that both spins 
components were included in (a), (b), and (e). 



If we compare with an aa ZGGNR, the absolute mag- 
netization increases to 1.73^s/cell and the total magne- 
tization to 1.0/is/cell in the a/3 ZGGNR. The side view 
of the spin density in Fig. @] shows that the /3-edge of the 
graphane nanoribbons is also spin polarized. 

The band structure of the a(3 ZGGNRs are shown in 
Fig. [5] The FM ground state shows a semiconducting be- 
havior [Fig.[B]Ja)] . There are two spin-polarized flat bands 
near the Fermi level with a direct bandgap of 0.58eV. 
Ferromagnetic semiconductors have attracted a great in- 
terest since 1990s because they may lead to spintronics 
applications, 2 - a/3 zigzag graphene/graphane nanorib- 
bons may be an alternative for the well-studied diluted 
magnetic semiconductors (DMS) materials^ with the 
fundamental difference that the system is purely two- 
dimensional in this case. The band structure correspond- 
ing to the nonmagnetic metastable state is represented in 
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FIG. 6: (Color online) The spin-resolved band structures of 
the aj3 zigzag graphene/graphane nanoribbons. (a) ferro- 
magnetic semiconducting ground state and (b) nonmagnetic 
metallic metastable state, (c) and (d) are zooms of (a) and (b) 
near the Fermi level, which was set to zero. The energy dif- 
ference between the ferromagnetic and the nonmagnetic state 
is Efm — Emm = 76meV/edge atom. 



Fig. lib). 

If we compare the energy difference between the mag- 
netic and the nonmagnetic cases for the three different 
types of nanoribbons shown in Table [IIJ we can see that 
the largest energy difference between the nonmagnetic 
and magnetic cases occurs for a/3 ZGGNRs. This indi- 
cates that the FM state of af3 ZGGNRs could be even 
more stable than the other cases considered in this work 
and summarized in Table UJJ 

a/3 ZGGNRs have an extra carbon atom and, there- 
fore, an odd number of atoms in the graphene nanorib- 
bons (13 carbon atoms in this case). This adds up one 
extra band near the Fermi level, which can be appreci- 



TABLE II: Electronic and magnetic properties of 
graphene/graphane nanoribbons with different hydro- 
genation at the interfaces. Bandgaps are in millielectron 
volt, energy differences (AE) are in millielectron volt per 
edge atom (meV/e.a.), and the total and absolute magnetic 
moments (/it and fi a bs, respectively) are in units of fiB /cell. 
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"We could not stabilize the AFM state for aB ZGGNRs. This 
could indicate that this state is not metaestable. 
6 The same situation ocurred for the FM state of 88 ZGGNRs. 



ated by comparing the nonmagnetic states of Figs. (He) 
and EJb) . The two bands that were previously dominat- 
ing the behavior near the Fermi level in aa ZGGNRs 
are now, respectively, pushed up and down in energy and 
very flat bands appear near the Fermi level, both for the 
FM ground state of Fig. UJa) and for the nonmagnetic 
metastable state of Fig. [6jb). 

Figures HJc) and [BJd) are zooms of the states near 
the Fermi level for the FM and non magnetic states, re- 
spectively. Even when the bands are very flat, there is 
a small dispersion that is caused by the interaction with 
the graphene/graphane interfaces. 

The results of the simple tight-binding model shown 
in Figs.|9jb) and[9je) reproduce qualitatively the behav- 
ior of the nonmagnetic band structure. However, the 
band at the Fermi level is completely flat. In order to 
obtain dispersion in this band, we need to include an 
interaction with the graphane nanoribbon. This can be 
done by introducing an extra row of carbon atoms cou- 
pled with hydrogen as in graphane, which simulates the 
graphene/graphane interface. In doing so we find an ex- 
tra band with a dispersion comparable to those from our 
ab initio simulations. The band width does not appre- 
ciably depend on the width of the graphene nanoribbon. 

More complicated tight-binding models for the case of 
an a(3 ZGGNRs where introduced in Ref. Ho]. However, 
the obtained band width-S was four times larger than the 
one found with our ab initio results. 

A very interesting feature of the PCD of the FM 
ground state of Figs. [5jc) and EJd) is that only one of 
the carbon sublattices is contributing to the edge state. 
The same happens for the nonmagnetic metastable state 
(not shown here). This can be incorporated in a simple 
tight-binding model represented in Fig. IHlJe) . All the co- 
efficients corresponding to one of the carbon sublattices 
are zero, B=D=F=H=J=L=0 in Fig.^e). 
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FIG. 7: (Color online) Side and perspective views of the 3D 
graph of the spin density of the antiferromagnetic ground 
state of the (3(3 graphene/graphane nanoribbons. The spin- 
polarized regions extend to include both edges of the graphane 
nanoribbon. Regions with positive spin density are repre- 
sented in red (dark gray) and regions with a negative spin 
density are shown in light blue (light gray). Isosurfaces cor- 
responds to (p-f — p^)=0.02 e/A 3 . 



Notice the nearly parabolic feature of the band of the 
nonmagnetic state of Fig.H^d) at k=0 and A;=0.5. This al- 
lows to define and effective mass (m e //) for the electrons, 
in this case m e ff — 5.77m e for k—0 and m e ff — 3.22m e 
for k=0.5, where m e is electron rest mass. 



C. Superlattice with the two edges 
half-hydrogenated: (3(3 zigzag graphene/graphane 
nanoribbons 
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FIG. 8: (Color online) The spin-resolved band structures of a 
(3(3 ZGGNRs. (a) Antiferromagnetic semiconducting ground 
state and (b) nonmagnetic semiconducting metastable state. 
The Fermi energy is set to zero. The energy difference be- 
tween the antiferromagnetic and the nonmagnetic state is 
Eafm — Enm = 37meV/edge atom. The spin-polarized 
bands coincide in both cases. 



We consider now the case of the (3(3 superlattice rep- 
resented in Fig. [He) ■ The spin-polarized ground state is 
the AFM state as was also the case in the aa superlattice. 
Figure [7] shows the spin density for the AFM ground state 
of the (3(3 ZGGNR. Both edges of the graphane nanorib- 
bon contributes to the spin polarization of the system. 
This can be better appreciated in the side view at the 
top of Fig. [7] The absolute magnetic moment of the 
AFM state is 1.92 /its/cell, almost twice larger than in 
the AFM state of the aa ZGGNR (1.09 /i B /cell). 

Figure [8ja) shows the band structure of the AFM 
ground state of the (3(3 ZGGNR and the nonmagnetic 
metastable state is shown in Fig. EJb). Both are semi- 
conducting states while the energy band gap of the AFM 
state (0.76eV) is larger than in the nonmagnetic case 
(0.29eV). It is interesting to note that even when we have 
a nonmagnetic state with a bandgap there is a transition 
to a magnetic state by further increasing the bandgap. 
The properties of both states are summarized in Table HT1 

The partial charge density graph of the AFM phase 
near the Fermi level obtained from Ef to Ef — 0.16eV 
shows that the edge states of the (3(3 superlattice extend 
also to the edge of the graphane nanoribbon as shown 
in Fig. [lje). The edge states decay more rapidly in the 
graphane nanoribbon than in the graphene region when 



we move away from the interface. 

The qualitative features of the (3(3 ZGGNRs are also 
captured by the tight-binding model of the Appendix as 
can be appreciated from Figs. HJc) andHJf). Figure [9j^f ) 
also correlates with the ab initio PCD of Fig.^c). This is 
discussed in more detail in the Appendix. We see that the 
different carbon sublattices contribute differently near 
each edge. This situation is similar to the case we ex- 
plained before for the aa ZGGNRs. 



III. CONCLUSIONS 

The degree of hydrogenation of the graphene /graphane 
interface has an important influence on the electronic 
and magnetic properties of the zigzag graphene/graphane 
nanoribbons. This parameter could be controlled exper- 
imentally to tune the bandgap and the magnetic proper- 
ties of graphene/graphane super lattices. 

We found that independent of the degree of hydrogena- 
tion of both interfaces the ground state of the system is 
semiconducting. However, the band structure near the 
Fermi level changes dramatically due to the contribution 
of the edge states in the three different super lattices. 
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While the ground state of the system is antiferromag- 
netic in aa and /3/3 ZGGNRs, it is ferromagnetic for a/3 
ZGGNRs. 

Our results show that aa ZGGNRs are a promising 
way to obtain stable graphene nanoribbons with similar 
electronic and magnetic properties as freestanding zigzag 
nanoribbons. The basic properties of freestanding zigzag 
graphene edges are preserved when the system binds on 
both edges to insulating graphane regions. We also found 
that the energy difference between zigzag and an arm- 
chair graphene/graphane nanoribbons is very small, the 
energy difference per atom is one order of magnitude 
smaller than the thermal energy at room temperature. 
This means that the formation energies of zigzag and 
armchair graphene/graphane nanoribbons are equally fa- 
vorable at room temperature. This situation is differ- 
ent from the case of freestanding nanoribbons where it 
was found that the zigzag edge is metastable, because 
it shows a planar reconstruction to lower its energy^ 
Then, graphene/graphane superlattices seem a promis- 
ing way to solve structural and stability problems faced 
by freestanding nanoribbons. 

We proposed a single atom thick ferromagnetic semi- 
conductor in a/3 zigzag graphene/graphane nanoribbons. 
It may be an alternative for the well-studied DMS 
materials^ 7 with the fundamental difference that the sys- 
tem is purely two-dimensional in this case. Diluted semi- 
conductors are of interest from both fundamental and 
applied points of view, and may soon lead to spintronics 
devices with new functionalities^ 

a/3 ZGGNRs can also be analyzed in a more general 
framework. In order to render a ZGNR half metallic, 
the symmetry between the left and right edges of the 
zigzag nanoribbon should be broken. For example, in 
Ref. [Tl| this was done by applying an electric field to 
the system, in Ref. |9| this was done by adding a Klein's 
edge to one side of the nanoribbon, and in Ref. [13 the 
two zigzag edges were terminated with different terminal 
groups, such as NO2 and CH3. The present work shows 
that in graphene/graphane superlattices, this could be 
achieved by controling the degree of hydrogenation at 
the two edges. 

From Fig. HJc) we can appreciate that the ferromag- 
netic ground state of a/3 ZGGNRs is almost a half metal. 
As the Fermi level is in between two spin-polarized bands, 
with little charge doping the Fermi level can be shifted 
either to the spin-up or spin-down band. This could be 
used in a potentially interesting electronic device. So 
half-metalicity in a/3 ZGGNRs can be achieved not with 
the help of an external electric held as in Ref. [ll| but 
with charge doping, which could offer new technological 
alternatives. 

It is known that the width of the nanoribbons has 
a strong influence on the bandgap of the graphene 
nanoribbons&ii as well as on the bandgap of the 
graphane nanoribbonsi 17 ' 18 Both, the width of the 
graphene nanoribbons and the level of hydrogenation at 
the interface can be used to control the electronic and 



magnetic properties of graphene/graphane nanoribbons. 
While the width of the graphane region can be used to 
control the energy window where graphene and graphane 
bands entangle. 

Most of the qualitative features obtained in our ab 
initio calculations can be explained with a simple tight- 
binding model incorporated in the appendix. This shows 
that the present results can be qualitatively extrapolated 
to other type of superlattices where graphene nanorib- 
bons are patterned by selective absorption of different 
atoms. Another interesting candidate, with similar or 
even better characteristics as graphane is graphene fluo- 
ride, a two-dimensional material that has been recently 
synthesized ! 22 ' 23 
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Appendix A: Tight-binding model of graphene 
nanoribbons 

In order to obtain extra insight in the pre- 
sented ab initio results for the electronic structure of 
graphene/graphane nanoribbons, we construct a simple 
tight-binding model. Here we focus only on the graphene 
part of the nanoribbon. To construct the tight-binding 
Hamiltonian we consider a tight-binding Bloch function 
per atom in the primitive cell 

I N 

v i 

where r*j is the position of a carbon atom in the primitive 
cell, and tp denotes the p z -like atomic wavefunction of a 
carbon atom. N is the number of considered unit cells. 
All carbon atoms are positioned on the (unrelaxed) hon- 
eycomb lattice as in pure graphene. We consider a zigzag 
nanoribbon with the same number N c of carbon atoms 
in the unit cell as in the graphene/graphane nanoribbons 
studied in this paper. For two a edges N c = 12 , for a 
nanoribbon with one a and one (3 edge N c — 13, and 
N c = 14 in case of two /3 edges. The total eigenfunction 
is then given by 

i=l 

This results in symmetric N c x N c tight-binding Hamil- 
tonians with the following non-zero off-diagonal elements 
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(with only nearest neighbor interactions) (1) nanoribbon 
with two a edges 

H(2i + l,2i + 2) = 7 cos(fcc) with i = . . . 5, (A3) 
H(2i + 2,2i + 3) = 7 with i = . . .4, (A4) 

(2) nanoribbon with one a and one f3 edge 

H(2i + l,2i + 2) = 7Cos(fcc) with i = . . . 5, (A5) 
H(2i + 2,2i + 3) = 7 with i = . . .5, (A6) 

(3) nanoribbon with two /J edges 

H"(2i + 1, 2i + 2) = 7 with £ = 0... 6, (A7) 
H(2i + 2, 2i + 3) = 7Cos(fcc) with £ = . . . 5.(A8) 



We have chosen the on-site energy as the zero of en- 
ergy J ip* (r — f*i)Htp(r — ri)dr — 0, and the tight-binding 
parameter 7 is defined in the usual way 



7 = / <p*{?— ?i)H(p(r- 



R)dr = 3eV (A9) 



where R is the vector connecting an atom with its nearest 
neighbor. 




H 1 1 1 1 1 \ 
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FIG. 9: (Color online) Band structures for graphene nanorib- 
bons with (a) two a edges containing 12 carbon atoms in the 
unit cell, (b) one a and one /3 edges and 13 atoms in the unit 
cell, and (c) two ft edges with 14 atoms in the unit cell. The 
Fermi energy is set to zero, (d), (e), and (f) show the val- 
ues of the coefficients \Ci\ 2 , for all atoms in the unit cell, of 
the highest valence band corresponding to (a), (b), and (c), 
respectively. 
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The band structures of the zigzag graphene nanorib- 
bons obtained within this tight-binding model are shown 
in Fig. |U Figures a)-[9fc) show results for the aa, a(3 
and /3/3 cases, respectively. The Fermi energy lies at 
Ep = 0. The obtained spectra agree qualitatively with 
the corresponding ab initio ones for the non-magnetic 
graphene/graphane nanoribbons in Figs. Etc), EJb), and 
IDJb). Note the absence of dispersion of the band at E = 
for the a/3 graphene nanoribbon [Fig. |9^b)] , in contrast 
to the ab initio result where a small dispersion was ob- 
served, see Fig. (6]Jd). This dispersion is a consequence 
of the graphane edge, which is not taken into account in 
the present tight-binding model. 

Below each band structure of Fig. [9l we plot the corre- 
sponding coefficients \Ci\ 2 from Eq. (|A2j) of the highest 
valence band eigenstate. Each curve is labeled by a letter 
indicating the corresponding atom in the nanoribbon unit 
cell. Figure EJd) shows that in case of an aa nanoribbon, 
the highest valence band state is located at both edges of 
the nanoribbon, represented by the carbon atoms A and 
L. The next most important contribution is not located 
at carbon atoms B and K, but at atoms C and J. Note 
that these two atoms (C and J) belong to different sub- 



lattices, and that C and A (and also J and L) belongs to 
the same sublattice. The results of Fig. [9jd) can be cor- 
related with the ab initio results of the PCD represented 
in Figs.[5ja) and[5jb). 

On the other hand, for the case of an a/3 nanoribbon 
[Fig. EJb)], the dispersionless band at the Fermi level is 
entirely built up from atoms that belong to one sublat- 
tice. The contributions coming from the other sublattice 
given by the coefficients B, D, F, H, J, and L are all zero. 
This is also the case for the ferromagnetic ground state 
of the aj3 nanoribbon as can be seen in the integrated 
local density of states of Figs. [SJc) andJSJd). 

The case of a nanoribbon with two /3 edges is similar to 
the aa nanoribbon. The highest valence band is again 
built up from the edge atoms M and N, see Fig. HHf). 
While the contribution from the neighboring atoms be- 
longing to the other sublattice (atoms A and L, or C and 
J) is again rather small. 

The present tight-binding model can be extended to 
describe the magnetic properties of the system, for ex- 
ample, by using the Hubbard model, as it was discussed 
in Ref . U 
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